Structural alphabets derived from attractors in conformational space

Background: The hierarchical and partially redundant nature of protein structures justifies the definition of frequently occurring conformations of short fragments as 'states'. Collections of selected representatives for these states define Structural Alphabets, describing the most typical local conformations within protein structures. These alphabets form a bridge between the string-oriented methods of sequence analysis and the coordinate-oriented methods of protein structure analysis.Results: A Structural Alphabet has been derived by clustering all four-residue fragments of a high-resolution subset of the protein data bank and extracting the high-density states as representative conformational states. Each fragment is uniquely defined by a set of three independent angles corresponding to its degrees of freedom, capturing in simple and intuitive terms the properties of the conformational space. The fragments of the Structural Alphabet are equivalent to the conformational attractors and therefore yield a most informative encoding of proteins. Proteins can be reconstructed within the experimental uncertainty in structure determination and ensembles of structures can be encoded with accuracy and robustness.Conclusions: The density-based Structural Alphabet provides a novel tool to describe local conformations and it is specifically suitable for application in studies of protein dynamics. © 2010 Pandini et al; licensee BioMed Central Ltd

GSATools: Analysis of allosteric communication and functional local motions using a structural alphabet

Motivation: GSATools is a free software package to analyze conformational ensembles and to detect functional motions in proteins by means of a structural alphabet. The software integrates with the widely used GROMACS simulation package and can generate a range of graphical outputs. Three applications can be supported: (i) investigation of the conformational variability of local structures; (ii) detection of allosteric communication; and (iii) identification of local regions that are critical for global functional motions. These analyses provide insights into the dynamics of proteins and allow for targeted design of functional mutants in theoretical and experimental studies. Availability: The C source code of the GSATools, along with a set of pre-compiled binaries, is freely available under GNU General Public License from http://mathbio.nimr.mrc.ac.uk/wiki/GSATools. Contact: alessandro.pandini@kcl. ac.uk or [email protected] Supplementary information: Supplementary data are available at Bioinformatics online. © 2013 The Author 2013. Published by Oxford University Press

Using local states to drive the sampling of global conformations in proteins

Conformational changes associated with protein function often occur beyond the time scale currently accessible to unbiased molecular dynamics (MD) simulations, so that different approaches have been developed to accelerate their sampling. Here we investigate how the knowledge of backbone conformations preferentially adopted by protein fragments, as contained in precalculated libraries known as structural alphabets (SA), can be used to explore the landscape of protein conformations in MD simulations. We find that (a) enhancing the sampling of native local states in both metadynamics and steered MD simulations allows the recovery of global folded states in small proteins; (b) folded states can still be recovered when the amount of information on the native local states is reduced by using a low-resolution version of the SA, where states are clustered into macrostates; and (c) sequences of SA states derived from collections of structural motifs can be used to sample alternative conformations of preselected protein regions. The present findings have potential impact on several applications, ranging from protein model refinement to protein folding and design.The computer resources, technical expertise, Red Española de Supercomputación, Cristian Micheletti, Michael Sadowski and Katarzyna Maksimia

Heart Failure Drug Modifies the Intrinsic Dynamics of the Pre-Power Stroke State of Cardiac Myosin

Omecamtiv mecarbil (OM), currently investigated for the treatment of heart failure, is the first example of a new class of drugs (cardiac myotropes) that can modify muscle contractility by directly targeting sarcomeric proteins. Using atomistic molecular dynamics simulations, we show that the binding of OM to the pre-power stroke state of cardiac myosin inhibits the functional motions of the protein and potentially affects Pi release from the nucleotide binding site. We also show that the changes in myosin ATPase activity induced by a set of OM analogues can be predicted from their relative affinity to the pre-power stroke state compared to the near rigor one, indicating that conformational selectivity plays an important role in determining the activity of these compounds

Molecular dynamics of electrosprayed water nanodroplets containing sodium bis(2-ethylhexyl)sulfosuccinate

The behavior of aqueous solutions of sodium bis(2-ethylhexyl)sulfosuccinate (AOTNa) subject to electrospray ionization (ESI) has been investigated by molecular dynamics (MD) simulations at three temperatures (350, 500 and 800 K). We consider several types of water nanodroplets containing AOTNa molecules and composed of a fixed number of water molecules (1000), N0 AOT AOT anions (N0 AOT = 0, 5, 10) and N0 Na sodium ions (N0 Na = 0, 5, 10, 15, 20): in a short time scale (less than 1 ns), the AOTNa molecules, initially forming direct micelles in the interior of the water nanodroplets, are observed in all cases to diffuse nearby the nanodroplet surface, so that the hydrophilic heads and sodium ions become surrounded by water molecules, whereas the alkyl chains lay at the droplet surface. Meanwhile, evaporation of water molecules and of solvated sodium ions occurs, leading to a decrease of the droplet size and charge. At 350 K, no ejection of neutral or charged surfactant molecules is observed, whereas at 500 K, some fragmentation occurs, and at 800 K, this event becomes more frequent. The interplay of all these processes, which depend on the values of temperature, N0 AOT and N0 Na eventually leads to anhydrous charged surfactant aggregates with prevalence of monocharged ones, in agreement with experimental results of ESI mass spectrometry. The quantitative analysis of theMD trajectories allows to evidencemolecular details potentially useful in designing future ESI experimental conditions

Binding Pockets in Proteins Induced by Mechanical Stress

We report the first observation of a pocket that opens as a result of a mechanical force applied to an Ig-like domain from the cardiac muscle. This previously unseen mechanism of pocket formation is revealed by molecular dynamics simulations under force. Preliminary investigations show that this “mechano-pocket” is potentially druggable and could be found in other domains from the same fold family, suggesting the existence of a general mechanism of pocket formation under mechanical stress

Reconstruction of ARNT PAS-B Unfolding Pathways by Steered Molecular Dynamics and Artificial Neural Networks

© 2021 The Authors. Several experimental studies indicated that large conformational changes, including partial domain unfolding, have a role in the functional mechanisms of the basic helix loop helix Per/ARNT/SIM (bHLH-PAS) transcription factors. Recently, single-molecule atomic force microscopy (AFM) revealed two distinct pathways for the mechanical unfolding of the ARNT PAS-B. In this work we used steered molecular dynamics simulations to gain new insights into this process at an atomistic level. To reconstruct and classify pathways sampled in multiple simulations, we designed an original approach based on the use of self-organizing maps (SOMs). This led us to identify two types of unfolding pathways for the ARNT PAS-B, which are in good agreement with the AFM findings. Analysis of average forces mapped on the SOM revealed a stable conformation of the PAS-B along one pathway, which represents a possible structural model for the intermediate state detected by AFM. The approach here proposed will facilitate the study of other signal transmission mechanisms involving the folding/unfolding of PAS domains.National Institutes of Environmental Health Sciences (R01-ES007685); Leverhulme Trust (RPG-2017-222)

Remote Working and Home Learning: How the Italian Academic Population Dealt with Changes Due to the COVID-19 Pandemic Lockdown

The COVID-19 pandemic introduced changes in people's lives that affected their mental health. Our study aimed to explore the level of psychological distress in the academic population during the lockdown period and investigate its association with the new working or studying conditions. The study sample included 9364 students and 2159 employees from five Italian universities from the study IO CONTO 2020. We applied linear regression models to investigate the association between home learning or remote working conditions and psychological distress, separately for students and employees. Psychological distress was assessed using the Hospital Anxiety and Depression Scale (HADS). In both students and employees, higher levels of distress were significantly associated with study/work-family conflicts, concerns about their future careers, and inadequacy of equipment; in employees, higher levels of distress were significantly associated with a lack of clarity on work objectives. Our results are in line with previous research on the impact of spaces and equipment in remote working/studying from home. Moreover, the study contributes to deepening the association between well-being and telework-family conflict, which in the literature is still equivocal. Practical implications require academic governance to promote sustainable environments both in remote and hybrid work conditions, by referring to a specific management by objectives approach

Specialized dynamical properties of promiscuous residues revealed by simulated conformational ensembles

The ability to interact with different partners is one of the most important features in proteins. Proteins that bind a large number of partners (hubs) have been often associated with intrinsic disorder. However, many examples exist of hubs with an ordered structure, and evidence of a general mechanism promoting promiscuity in ordered proteins is still elusive. An intriguing hypothesis is that promiscuous binding sites have specific dynamical properties, distinct from the rest of the interface and pre-existing in the protein isolated state. Here, we present the first comprehensive study of the intrinsic dynamics of promiscuous residues in a large protein data set. Different computational methods, from coarse-grained elastic models to geometry-based sampling methods and to full-atom Molecular Dynamics simulations, were used to generate conformational ensembles for the isolated proteins. The flexibility and dynamic correlations of interface residues with a different degree of binding promiscuity were calculated and compared considering side chain and backbone motions, the latter both on a local and on a global scale. The study revealed that (a) promiscuous residues tend to be more flexible than nonpromiscuous ones, (b) this additional flexibility has a higher degree of organization, and (c) evolutionary conservation and binding promiscuity have opposite effects on intrinsic dynamics. Findings on simulated ensembles were also validated on ensembles of experimental structures extracted from the Protein Data Bank (PDB). Additionally, the low occurrence of single nucleotide polymorphisms observed for promiscuous residues indicated a tendency to preserve binding diversity at these positions. A case study on two ubiquitin-like proteins exemplifies how binding promiscuity in evolutionary related proteins can be modulated by the fine-tuning of the interface dynamics. The interplay between promiscuity and flexibility highlighted here can inspire new directions in protein-protein interaction prediction and design methods. © 2013 American Chemical Society